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Abstract 

The Standard Solar Model (BP04) predicts a total neutrino flux that is 17.2% larger than 
measured in the salt phase of the SNO detector (and if it were significant it will indicate oscillation 
to sterile neutrinos). Hence it is important to examine in details uncertainties (and values) of 
inputs to the SSM. Currently, the largest fractional uncertainty is due to the new evaluation of 
the surface composition of the sun. We examine the nuclear input on the formation of solar 
[5*17(0)] and demonstrate that it is still quite uncertain due to ill known slope of the measured 
astrophysical cross section factor and thus ill defined extrapolation to zero energy. This yields 
an additional reasonably estimated uncertainty due to extrapolation of I3Q eV-b {^i^yj- Since a 
large discrepancy exists among measured as well as among predicted slopes, the value of 6*17(0) is 
dependent on the choice of data and theory used to extrapolate 5*17(0). This situation must be 
alleviated by new measurement (s). The "world average" is driven by the Seattle result due to the 
very small quoted uncertainty, which we however demonstrate it to be an over-estimated accuracy. 
We propose more realistic error bars for the Seattle results based on the published Seattle data. 



PACS numbers: 25.20.Dc, 25.70.De, 95.30-K, 26.30.+K, 26.65.+t 
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I. INTRODUCTION 



The high precision measurement of neutral current interactions of "^B solar neutrinos" in 
the SNO detector (with added salt) yields the measured flux: $7vc = 4.94 ± 0.21 (stat) 
-a34 (syst) X 10^ cm~^ sec"-^, with a total uncertainty of +8.8 -8.1%. The most realistic 
Standard Solar Model (SSM) prediction for the flux, labeled as BP (04) in is 5.79 x 
10^ cm~^ sec~^, hence: 

|^(z/) - 1 = 17.2% (1) 

It is important to examine if this discrepancy is significant as it could for example indicate 
(further) oscillation of (solar) electron neutrinos to sterile neutrinos. The uncertainties of 
the SSM prediction of solar neutrino fluxes were studied by Bahcall and Serenelh as listed 
in Table 1. 

TABLE I: Fractional uncertainties of predicted and ^Be solar neutrino fluxes 



Source 


^B 




PP{Sn) 


1% 


0.4% 


^He +3 He{S^s) 


2% 


2% 




8% 


8% 


Be{Sn) 


4% 





Composition (Z/X) 


12% 


5% 


Opacity (k) 


5% 


3% 


Diffusion 


4% 


2% 


Luminosity {Lq) 


3% 


1% 


Total: 


16.7% 


10.4% 



Currently the largest fractional uncertainty of the ^B solar neutrino flux is due to the 
chemical composition CZi/X) which is a consequence of the new evaluation of the abundance 
of heavy elements (mainly oxygen and neon) on the surface of the sun. The new evaluation 
leads to doubling the uncertainty (from 6% to 12%) and to what was recognized as a new 
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"standard solar model problem" vis-a-vis the predicted convective ^^^^ ^> shown 
to be in conflict with SSM predictions based on helioseismology 5]. Hopefully this large 
uncertainty due to Z/X will be resolved soon, and there is progress in reducing the error 
due to 5*34(0), see Table, leading to a considerably smaller total error. But in the same 
time the errors of nuclear inputs are far from resolved. In this paper we address the error 
of a crucial nuclear input, the formation of solar via the '^Be{p, j)^B reaction and its 
associated astrophysical cross section factor at zero energy, 5*17 (0). 

The value of 5*17(0) = 21.4 eV-b and the (impressive) accuracy of ±4% listed in Table 1 is 
adopted from the Seattle^roup 0| . It should be compared to the previously quoted value of 
5*17(0) = 19 +4 -2 eV-b ^7,]- The new quoted value of 5*17(0) leads to an increase of 12.6% in 
the predicted ^B solar neutrino flux. We examine in details the claims of accuracy as well as 
the quoted value. We demonstrate that the "world average" value of 5*17(0) is driven by the 
Seattle result with an over estimated accuracy, and more realistic error bars are proposed. 
Due to large discrepancy between measured slopes of (mostly DC) data, and since theories 
also disagree on the predicted slopes, the value of 5*17(0) is dependent on the choice of data 
and theory that one uses to extrapolate 5*17(0). This situation must be alleviated by new 
measurement (s) . 



A. Direct Capture 



It is well known that the "old" data on Direct Capture (DC) Be{p,'j)^B reaction 
y, 10,^1 exhibit major systematic disagreements. But the situation is not improved with 

Q|, as we show in Fig. 1. The data of the Orsay group 



"modern" data on DC |6| , |12l. 113 

jl^ and Bochum group do not agree with that of the Seattle group [6] and Weizmann 
group ^|. As can be seen from Fig. 1 the disagreements among individual 5*17 data points 
are by as much as five sigma and there is not a single measured data point of the Bochum 
group that agrees with a data point measured by the Seattle group. Most disturbing is the 
systematic disagreement at low energies of the two "modern" high precision measurements 
of the Seattle p and Weizmann groups. These two measurement are in good agreement 
at higher energies (around 1 MeV), but the Weizmann data are systematically below the 
Seattle data at low energies (300 - 500 keV), indicating a different slope. 

The large discrepancies between measured individual DC data points shown in Fig. 1, 
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FIG. 1: Modern Direct Capture world data on Sn (excluding data around the 632 keV resonance) 
exhibiting marked systematical deviations by up to five sigma. 

demonstrate large systematical differences between "modern" DC data. These large sys- 
tematical differences must be resolved before these data are used to extract 5*17 (0) from DC 
world data. The systematically disagreeing DC data can not be handled algebraically using 
statistical methods to extract average properties of DC data and in this case chi-square tests 
are meaningless. 



II. THE SLOPE MEASURED BETWEEN 300 - 1,400 KEV 



The slope of the measured astrophysical cross section factor between 300 and 1400 keV is 
essential for extracting the d-wave contribution to measured data. This contribution must 
be subtracted from the data in order to extrapolate to zero energy ^ where the s-wave 
dominate (approximately 95%). But so far we can only rely on theoretical estimates of 



the d-wave component |l5l . |l 
measurement of the slope. 



mm 

and these models can be tested only by a precise 
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FIG. 2: The slopes (S' = dS/dE) of world data measured between 300 and 1400 keV. All available 
DC data are shown and the slope of the RIKEN2 data 2] is shown correctly. The shown range 
of so called "average values" for DC and CD are in agreement as discussed in the text. 

The Seattle group extracted the average slope of DC data with high pree,s,o,r (±4.5%), 
as shown in Fig. 19 of their paper pj. But three of the five DC data points selected to be 
shown in Fig. 19 ^6] (note that not all available DC data are shown in Fig. 19) have central 
values that disagree by a factor of 2.5 and almost three sigma. Furthermore, the quoted 
high precision is in fact driven by the claimed high precision of Seattle experiment hence in 
the next section we examine this precision in details. 

In Fig. 2 we show the slope of all available DC data measured between 300 - 1400 keV, 



including the data of Vaughn Parker jicl], and Kavana gh |ll| tha t were omitted in Fig. 

|21|, 



19 of [6]. We also show all thus far measured CD data [l^l, |2m21|, |2^ . We conclude that the 
slope parameter can not be extracted from DC data with the (impressive) accuracy of 4.5% 
0], unless one excludes some of the DC measurements discussed above. An error which is at 
least a factor of 2 larger seems like a more reasonable choice. In the same time it is obvious 
from Fig. 2 that there is no evidence for large disagreement of the slopes of CD and DC 
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data as claimed in p]. In any case the large dispersion of measured slopes shown in Fig. 2 
can not allow us to determine the slope from either CD or DC data with sufficient precision 
and this question must be addressed by future experiments. Lack of accurate determination 
of the slope must lead to an additional uncertainty due to the ill defined subtraction of the 
d-wave and the extrapolation to zero energy. 



III. ACCURACY OF THE SEATTLE DATA 



Since the claimed high precision of the Seattle measured data points and quoted extracted 
5*17(0) 0] drive the value and the error of so called "world average" it seems reasonable to 
examine it in detail, especially at low energies. In Fig. 3 we show the results of the target- 
beam calibration data shown by the Seattle group using the ''i?e(a, 7)^^C reaction. As can 
be seen in Fig. 3 both calibration spectra published in the Phys. Rev. Lett, and the Journal 
Phys. Rev. C (which where measured with different Be targets) exhibit a resonance energy 
which is off by 9 keV from the known energy in ^^C . This 9 keV shift is intolerable and 
the authors as yet did not explain this discrepancy in a published erratum. The repeated 
publication of the same mistake first in the PRL paper and a year later in the PRC long 
paper does not lend credence to the claimed high precision experiment. 

Furthermore, in their paper they show the ^ Be target profile measured in the middle 
(PF2) and at the end (PF3) of the experiment j^. These profiles are superimposed on top 
of each other in Fig. 4a, from which it becomes clear that we can not support the claim 6| 
that "they are similar". In fact there is evidence that the Be moved further into the target 
in the intervening period between these two measurements. The unstable nature of the ^ Be 
during the experiment must be refiected in a systematic uncertainty due to variation of the 
expected yield, as we discuss below. Another systematic uncertainty that is not considered 
here is associated with the evaluation of the effective center of mass energy. 

These profiles are measured with alpha-particle beams with energies of approxi- 
mately 1.4 MeV in steps of AE = - E/j which can be translated to steps of AE 
= (E„ - E^) X #p+(M„+7Be)/#a+(Mo+TBe) P^^^^^ ^cams. Note that the ratio 
^p+(Mrr^'^Be)l ^a+{Mo+'^Be) ^Imost independent of the exact Mo + ^ Be mixture of the 
target jy|, as well as variation among the various tabulated individual energy loss. Hence 
for an incoming proton beam (Ep) the Be target profiles during the Be{p,'y)^B reaction 
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FIG. 3: The measured beam-target calibration spectra a| using the ^i?e(a,7)^^C reaction. Both 
spectra pubhshed in PRL and PRC exhibit resonance energy which is off by 9 keV from the well 
known resonance energy in ^^C as indicated. 



measurement is given by: 



dNt 



dE 



- {E^ - 1370) X J^2±iM2±l^) oc Profile(K 

dX a + (Mo+T B<i) 



(equ. 2) 



The so obtained profile is shown in Fig. 4b, it leads to very large target thickness, for 
example of the order of 35 keV at a beam energy of approximately 130 keV. Over such 
a thick target the cross section varies by almost a factor of five, as shown in Fig. 4c. A 
convolution of the so obtained target profile with the expected variation of the cross section 
leads to the expected yields: 



Yield oc xa(^- 



(equ. 3) 



The yields evaluated for the two shown profiles measured in the middle (PF2) and at the 
end (PF3) of the experiments differ by 7.5% at the lowest measured energy, almost a factor of 
3 larger than the quoted systematic uncertainty 0]. In Fig. 5 we show the differential yields 
= Y(PF2)/Y(PF3) - 1 evaluated for all incoming proton beam energies. This differential 
yield must be considered a lower limit of the systematical uncertainty. This uncertainty is 
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FIG. 4: (a) Target profiles measured using tlie Be{a,j) C reaction in the middle (PF2) and 



at the end (PF3) of the experiment 



, (b) the target profile for a 150 keV proton beam, (c) the 



predicted cross section variation across the target. 
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FIG. 5: One systematic error of the Seattle experiment ^ due to variation in ^Be target profile 
shown in Fig. 4. The two profiles measured in the middle (PF2) and at the end (PF3) of the 
experiment are used. 
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FIG. 7: A comparison of the Descouvemont cluster model VS\ and Typel potential model 
with the GSI 12^1 and Weizmann data Q- 

a factor of 3 larger than the quoted varying systematical uncertainty 0|. Clearly the most 
troubling uncertainty is at energies below 400 keV. This energy range was emphasized in 
the Seattle paper as most useful for extracting <S'i7(0), but it is measured with the largest 
systematical uncertainty. 

The large systematical uncertainty of Seattle data at low energies is particularly bother- 
some vis-a-vis the exclusion of theoretical models such as the potential model of Typel 2^ ■ 
In Fig. 6 we show a comparison of that model to the Seattle data where we find that the 
model reproduce well the high energy data but not the low energy data (or vice versa). The 
large systematical error at low energies shown in Fig. 5 cast doubt on the claim that the 
Typel model is "rejected" by the data. This model fits well all other data and consistently 
yield extrapolated values of S'i7(0) which are approximately 3 eV-b smaller, as shown in Fig. 
7. 

We conclude that in the absence of an accurate measurement of the slope we must add 
an extrapolation error to the quoted S'i7(0). And in view of the fact that the Typel model is 
consistent with all existing data except the Seattle data, and in view of the large systematical 
error of the Seattle data at low energies, a reasonable estimate of an extrapolation error is 
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of +0.0 -3.0 eV-b. Such an a large error implies the discrepancy of the SSM shown in equ. 
1 is not significant solely due to nuclear uncertainty even if the uncertainty due to chemical 
composition will be reduced. 
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